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Results are presented from the numerical calculations 
carried out to produce detailed information on the kinematic, 
dynamic, dissipative, and thermodynamic characteristics of a 
uniform half infinite stream , mixing with a quiescent fluid 
of the same composition. 

An effective Prandtl number of Pr^ = 1 is assumed for the" 
constant pressure, non-isoenergetic turbulent mixing process. 

The ratio of specific heat becomes absorbed by selecting 
Crocco number instead of Mach number as measure for the com- 
pressibility, thus allowing generalization of the results to 
any perfect non-reacting gas having constant specific heat. 

After identifying a functional form for the streamwise 
component of the velocity profile , a single empirical mixing 
parameter becomes well defined and can be absorbed in a rational 
presentation of structural details of jet mixing regions, such 
as the vertical velocity component, temperature and density dis- 
tribution, integrals describing flow of mass, momentum, mechan- 
ical energy, the transfer of shear work and heat across individ- 
ual streamlines, as well as local and integrated dissipation 
rates for mechanical energy. 

Information on the empirical parameter a remains generally 
incomplete. Although values for low-speed isoenergetic mixing 
are well established and effects of Mach number have been ten- 
tatively reported, no such information is presently available for 
temperature level and temperature differential influence. 
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Introduction 


The interest in flow problems involving separation from solid 
boundaries leads to studies of the flow mechanism in the wake, and its 
component flow regions. 

Significant wake flow components may be identified as: 

i) the flow field near the separation point, 

ii) the jet mixing component between the wake and the adjacent 
free stream, 

iii) the region of reattachment (often recompression) at the end 
of the wake, 

iv) the flow field within the wake, 

v) the redevelopment of the flow field downstream of the end of 
the wake . 

Understanding of wake dynamics has benefited greatly from the 
analysis of such wake flow components, as evidenced by the treatment of 
the classical base pressure problem and its ramifications (ref. l). 

The jet mixing component is of special importance in its contribution 
to wake flow mechanisms, and has, therefore, received wide attention in 
the literature . 

Mixing between a uniform stream and a quiescent fluid, as well as 
mixing involving two streams have been studied extensively both analyti- 
cally and experimentally.* Yet, even, the most fundamental and simplest 
information, e.g. concerning the influence of compressibility and 
temperature level on the effective turbulent eddy viscosity in turbulent 
jet mixing regions is still subject to searching speculations (ref. 4). 
While it appears logical to draw such information from an analysis of 
experimental data on mixing profiles, in concentrating on the vicinity 
of the inflection point tangent, and to observe the close relationships 
between the functional presentation of a mixing profile and the result- 
ing definition of an empirical mixing parameter, such an approach has, 
so far, not been uniformly accepted. After identifying a functional 
form for the velocity profile, such a single empirical parameter becomes 
well defined and can be utilized in a rational presentation even of de- 
tailed information on jet mixing regions. Final identification by 
numerical values may still have to depend on further empirical informa- 
tion on this mixing parameter. The present communication is concerned 
with the turbulent, compressible constant pressure mixing problem, 


*The rapid rate at which contributions are being made to these 
fields forces one to look beyond the treatments and references included 
in such standard works as by Schlichting (ref. 2) and Pai (ref. 3) and- 
to survey the current literature . 
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resulting from interaction between a uniform stream and a quiescent wake, 
both having the same composition, but, in general, different stagnation 
temperature and an effective turbulent Prandtl Number of unity. De- 
tailed information is given on the kinematic, dynamic, dissipative and 
thermodynamic structure of such jet mixing regions. The use of Crocco 
Number as parameter for compressibility effects eliminates the influence 
of the specific heat ratio, and thus, permits a wider utilization of 
the results (ref. 5)* It shall be noted that jet mixing between a 
uniform stream and a quiescent fluid serves satisfactorily as a flow com- 
ponent in many problems concerned with wake dynamics. Thermodynamic 
analysis of separated flow regions and a study of the mechanism of energy 
transfer to and across wakes requires however the consideration of finite 
entrainment velocities in the jet mixing component (ref. 6). A subsequent 
report similar in scope to the present one, will present information on 
the two-stream jet mixing problem.* 

In the study of the flow mechanism and heat transfer in separated 
flow problems, it is necessary to study the detailed local properties 
within such turbulent jet mixing regions which would be helpful to 
understand the basic flow mechanisms and identify the controlling com- 
ponents within such separated flow regions. 

The present work is essentially a logic extension of previous 
investigations conducted at the University of Illinois, which utilized 
an integral momentum approach to obtain theoretical solutions for the 
isoenergetic (ref. 5> 7) and non-isoenergetic (ref. 8, 9) compressible 
turbulent jet mixing problem. In order to have consistent presentation, 
the previous work on which the present analysis is based will be first 
summarized and reviewed. 


*A comprehensive computation program has been carried out at the 
University of Illinois under the NASA Grant NsG-13-59 and samples have 
been presented in ref. 6. 
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Theoretical Analysis 


The differential equation of motion for a constant pressure 
turbulent jet mixing region was highly simplified and solved through 
an integral transformation by extending Prandtl’s exchange coefficient 
concept. The solution was written as functions of the initial dis- 
turbed profile as well as the "position parameter Tp" . (Ref. 5). It 
was shown that in the case of small initial disturbance or at a location 
far downstream, the velocity profile would asymptotically reach one 
termed as "fully developed", (^p = 0) which is no longer dependent upon 
the initial disturbed profile. 

For a turbulent Prandtl Number of one (Pr t =1), the Crocco 
integral energy relationship will relate the stagnation temperature 
profile throughout such a mixing region uniquely to the velocity pro- 
file. (Ref. 9). 

The solution for the flow field in the mixing region thus obtained 
was interpreted to hold in an intrinsic coordinates system (x,y) which 
was subsequently localized with respect to the reference coordinate 
system (X, Y) by a momentum integral relationship. In the following, 
all equations are written for the "fully developed" profiles within 
such a constant pressure compressible non-isoenergetic turbulent jet 
mixing region. 


Jet Mixing Profile 


The dimensionless velocity profile is given by 


where 


o = \ (l + erf T]) 



T1 - a X 

X 


( 1 ) 


and ct is the similarity parameter for the homogeneous coordinate system. 
A discussion of the parameter a as identified for the error function 
velocity distribution is given in Appendix A. 
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( 2 ) 


The stagnation temperature profile is given by 
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The dimensionless shift of the "intrinsic system of coordin- 
ates" with respect to the reference coordinate system is given by 
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where ^ is a large value of T| such that 


1 - 9 (%) < t 

and | 1 - A (T^) | < t' 

t and t' being small quantities. 


Auxiliary Integrals 

The integrals related to this analysis are defined and listed as 
follows: 
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*Note that this integral 1^ has been defined as 1^ in Refs. 1 and 9* 
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For fixed values of C 2 and — , these integrals may be repres- 

a oa 

ented in short by I x (Tl), I Q (T|), I ? (T|) and I 4 (Tl) respectively. 


Jet Boundary Streamline 

The "jet boundary streamline” which separates the fluid of the 
external stream from the fluid entrained within the wake is identified 
by T|j which satisfies: 

II (%) - Ii (Tlj ) = I 2 OTr) ( 4 ) 


Energy Transfer 

The energy transferred across the "jet boundary streamline" is 
given by q which satisfies 

= St X a = (1 - C 2 ) E/l - ( 5 ) 

C p u (T - T, ) oa 

p a a oa d 

where 

E E [ U (c/ . 

- [ x 3 K > 

All the integrals and quantities mentioned above have been 
calculated on the Illiac* and have been tabulated and presented in 
graphical form (Ref. 7 an d 9 ) • 
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Local Characteristics within the Jet Mixing Region 

In view of the fact that the treatment presented above was based on 
integral relations, lateral differentiations of the profiles to obtain 


*Electronic digital computer, Engineering Research laboratory, 
University of Illinois. 
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local properties are not recommended. The present analysis is again 
based on integral relations whenever possible. 

Following a certain streamline at the ordinate y within such a 
jet mixing region, one would have from the continuity relationship that 
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After the differentiation is carried out, the equation (6) goes 
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which is a basic relation needed in this analysis. It describes the 
change of the dimensionless coordinate 71 pertaining to a certain stream 
filament as it proceeds downstream. It is obvious that the jet boun- 
dary streamline has the constant property values within the "fully de- 
veloped mixing regions. 


Selecting a control volume within such a jet ^nixing region such 
that the top control surface coincides with a streamline (see Fig. l), 
the integral momentum relationship 

y y x 
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becomes after being differentiated with respect to x 
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which is equivalent to 


T t = 


(1-C 2 )p u 2 , 
v a a a d 

a dx 


J 




_oo A _c 2 cp 2 
a T 


n 

dTl-x r 

J n 


_2_ 


t . A -C 2 tp 2 
J a 


dT> 


( 8 ) 


Carrying out the differentiation, the equation (8) goes into 
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Defining the shear stress function t by 
t = T t a 


(9) 


p U 

a a 

one obtains by substituting equation (T) into the equation (9): 

T = (1 - C 2 ) [I (Tl) - cp (i(-n) - I (T.))] 
a 2 1 1 j 

It is readily seen from the definition 

„ . du p a a dtp 

t >y p x dt 
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one can estimate the eddy diffusity e with the help of equation (10) 
and obtain 


(10) 


£ = = (A -C 2 <p 2 ) /tt [I (Tl) - cp (I (Tl) - I (Tl.))] 

a 2 1 1 J 


XU 


(ID 


which shows that at a location x, e is certainly not constant along the 
y direction. Limited experimental results obtained for C a 2 = 0 (Ref. 10) 
seem to support our analysis, see Figure 2a. While no experimental veri- 
fication of the Mach Number influence or the temperature ratio effect is 
presently available, it is reasonable to assume that our analysis should 
again exhibit the correct trends.* 


The derivative of the shear stress function is also given as 
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By use of the Reynold’s Analogy for the case of unity turbulent 
Prandtl Number , 


*It must be noted that our results are in strong disagreement 
with the speculations of Ting and Libby(Ref. 4). 
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A direct derivation of equation (13) by applying the energy 
balance is shown in the Appendix B. 

The local shear work and dissipation functions, W g , 9^ can be 
evaluated from the following expressions namely 
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which represents the energy transfer across the streams and has 
been presented before (equation 2-3, in Ref. l). 
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The integral 
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dissipation up to ordinate y in the mixing region and it is shown in 
the Appendix C that 
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For y approaches infinity, one would obtain i> , the total time rate of 
dissipation of mechanical energy per unit length along the mixing region 
(with unit thickness) which is given by 
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I (TU) = I (TL) - I (7 ) 
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has been introduced. 


The velocity component v in the y direction can also be obtained 
for such a jet mixing region. In Appendix D, it is shown that the v 
component measured in the reference system of coordinates (X, Y) is 
given by 


a cp> = cp (VHj - [I (Tl) - I ('ll,)] 


(18) 


where 


cp 1 = 


u 


*Note that this is essentially the integral mechanical energy 
relationship 



which can be derived by integrating along the y direction the equation 
of motion for such a constant pressure jet mixing region after being 
multiplied by the velocity. 
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Numerical Results 


Evaluations of all these functions and integrals introduced 
above have been carried out on the Illiac for five parametric values 

T b T b 

rji > (ip — 0*1^ 0. 1^ 5 ) 3-t five values of Crocco Number 

oa oa 

(c a 2 = 0, 0.2, 0.4, 0.6, 0.8). 

The distribution of the shear stress function t and the eddy 
diffusivity function e are presented in Fig. 2 for three values of 

— at various Crocco Numbers. In Fig. 2a, experimental results of 
oa 

T b 

limited amount for isoenergetic (=p- = 1) incompressible flow (C 2 = 0) 

ca a 

obtained by hot wire measurements (ref. 6) are also shown. The agree- 
ment between the theoretical calculations and the experimental results 
are reasonably good. 

The distribution of shear work and dissipations functions W 9 

s D 

are presented in Fig. 3. The total time rate of dissipation of mechan- 
ical energy in such jet mixing regions is presented in Fig. 4 while the 
distribution of vertical component of the velocity (v- component ) within 
the jet mixing region is presented in Fig. 5 . 
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Appendix A 

Discussion of the Similarity Parameter o 
in the Turbulent Jet Mixing Region 


In a fully-developed jet mixing region with Pr^ = 1, the Similarity 
parameter a can best be related to the change in slope of the inflection 
point tangent of the velocity profile. Selecting the established "erf- 
function" profile in such a mixing region, one can obtain by differenti- 
ation 


dcp 

dil 



e 


-7f 


so that for the inflection point of the velocity profile, 71 ■ 0, one 
would obtain 

*0 = i 

d V T\ = 0 ^ 


Within the x, y system of coordinates one may determine, for the 
stations, x^ and x 2 the maximum slope of the velocity profile, 
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As x^ becomes large, a stationary value is obtained for a (fully 
developed profile ) 

/tT 


u 


5u 

3y 


The value of o has been well established for the isoenergetic 
incompressible flow case to be a = 12, and the relation of a = 12 + 

2. 758m has been also suggested by Korst and Tripp (ref. 11) to account 
for the effect of compressibility. 
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Appendix B 

Energy Transport within the Turbulent Jet Mixing Region 



Selecting a control volume as shown in the Fig. B-l such that the top 
control surface coincides with a streamline considered, the integral 
energy relationship 
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which is exactly the Equ. (13). 
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Appendix C 

Derivation of Equation (l6) 


Starting from the Equ. (10) 
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Also 
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Equation (C-l) can be brought into the form 
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Appendix D 

The v-Coraponent of the Velocity within the Jet Mixing Region 


From the defining relationship of the intrinsic coordinate 
system, namely 


X = x 

Y = y - y m 


one obtains by following a certain fluid element in the reference 
system of coordinates 


As 
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one combines the Equations (D-l) (D-2) (D-3 ) (D-4) with the 
additional relationship given by Equation (7) and obtains 

a cp* = 9 (1) - \) - III (n) - lx (Tlj)] (A-C a V)* (D-5) 

The distributions of the v-component of the velocity are plotted 
in Figures 5a,b,c for various values of free stream Crocco Numbers 
and stagnation temperature ratios across the mixing region. 


* Note that the dimensionless velocity components for the x, y 
system and the X, Y system are related by 

cp(X, Y) = ep(x,y) 

cp*(X,Y) = cp*(x,y) - cp(x,y) 7^ 
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Fig. 5c The Distribution of the v-Component of the Velocity for T /T * 2.0 





